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Abstract. Designed as a high-sensitivity gamma- 
ray observatory, the Fermi Large Area Telescope is 
also an electron detector with a large acceptance 
exceeding 2 m 2 sr at 300 GeV. Building on the 
gamma-ray analysis, we have developed an efficient 
electron detection strategy which provides sufficient 
background rejection for measurement of the steeply- 
falling electron spectrum up to 1 TeV. Our high 
precision data show that the electron spectrum falls 
with energy as E ~ 3 and does not exhibit prominent 
spectral features. 
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I. Introduction 

Accurate measurements of high-energy Cosmic Ray 
(CR) electrons (we hereafter refer to electrons as a sum 
of e + and e~ unless specified otherwise) provide a 
unique opportunity to probe the origin and propagation 
of CRs in the local interstellar medium and constrain 
models of the diffuse gamma-ray emission [1]. Prior to 
2008, the high-energy electron spectrum was measured 
by balloon-borne experiments [2] and by a single space 
mission (AMS-01, [3]). The measured fluxes differ by 
factors of 2 — 3. While these data allowed significant 
steps forward in understanding CR origin and propaga- 
tion, constraints on current models remain weak. The CR 
propagation package GALPROP [4], on the assumption 
that electrons originate from a distribution of distant 
sources mainly associated with supernova remnants and 
pulsars, predicts a featureless spectrum from 10 GeV 
up to few hundreds of GeV. Above that energy, due 
to the actual stochastic nature of electron sources in 
space and time, and to the increasing synchrotron and 
inverse Compton energy losses, the spectral shape may 
exhibit spatial variations on a scale of a few hundred 
parsecs. Nearby sources start contributing significantly 
to the observed local flux and may induce important 
deviations from a simple power law spectrum [2], [5], 
[6]. 

Recently published results from Pamela [7], ATIC [8], 
H.E.S.S. [9] and PPB-BETS [10] have opened a new 
phase in the study of high energy CR electrons with a 
new generation of instruments. These four experiments 
report deviations from the reference model mentioned 
above. Pamela measures an increase of positrons with 
respect to electrons at energies above a few GeV; ATIC 
and PPB-BETS detect a prominent spectral feature at 
around 500 GeV in the total electron plus positron 



spectrum; H.E.S.S. reports significant steepening of the 
spectrum above 600 GeV. All these results may indicate 
the presence of a nearby primary source of electrons 
and positrons. The natures of possible sources have been 
widely discussed, two classes of which stand out: nearby 
pulsar(s) ([11], [12] and references therein) and dark 
matter annihilation in the Galactic halo, e.g. [13]. 

The Large Area Telescope (LAT) is the main instru- 
ment on-board the Fermi Gamma-Ray Space Telescope 
mission. It is conceived as a multi-purpose observatory 
to survey the variable gamma-ray sky between 20 MeV 
and 300 GeV, including the largely unexplored energy 
window above 10 GeV. It is therefore designed as a 
low aspect ratio, large area pair conversion telescope to 
maximize its field of view and effective area. The LAT 
angular, energy and timing resolution rely on modern 
solid state detectors and electronics [14]: a multi-layer 
silicon-strip tracker (TKR), interleaved with tungsten 
converters with a total depth of 1.5 radiation length (Xo) 
on-axis; a hodoscopic CsI(Tl) calorimeter (CAL), 8.6 
Xo deep on-axis; a segmented anticoincidence plastic 
scintillator detector (ACD) and a flexible, programmable 
trigger and filter logic for on-board event filtering. 

Since electromagnetic (EM) cascades are germane 
to both electron and photon interactions in matter, the 
LAT is also by its nature a detector for electrons and 
positrons. Its potential for making systematics-limited 
measurements of CR electrons was recognized during 
the initial phases of the LAT design [15], [16] and 
demonstrated by the recent publication of the Fermi 
measurement of the CR electron spectrum between 20 
GeV and 1 TeV [17]. In this paper we describe the 
technique developed for this measurement and our vali- 
dations using ground and flight data, and present the first 
high statistics CR electron spectrum from 20 GeV to 1 
TeV based on the data taken in the first six months of 
the mission. A more complete description of the analysis 
procedure is available in other contributions from the 
Fermi collaboration to these proceedings ([18], [19]). 

II. Event Selection 

An analysis that separates electrons from the dominant 
CR hadrons is required to fully exploit LAT's large 
collecting power and long observation time. The on- 
board filter is configured to accept all events that deposit 
at least 20 GeV in the calorimeter; thus we ensure 
that the rare high-energy events, including electrons, 
are available for thorough analysis on the ground. We 
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developed a dedicated event selection for high energy 
electrons that provides a large geometry factor with 
a residual hadron contamination less than 20% at the 
highest energy. At any given energy, the geometry factor 
(GF) is defined as the proportionality factor relating the 
rate of events passing the selection criteria to the incident 
flux. It measures the instrument acceptance, i.e the 
integral of the effective area over the instrument field of 
view. It is numerically evaluated using a MC data sample 
of pure electrons. As for the analysis developed for 
extracting LAT photon data [14], the electron selection 
essentially relies on the LAT capability to discriminate 
EM and hadronic showers based on their longitudinal 
and lateral development, as measured by both the TKR 
and CAL detectors. The background rejection power for 
photon science is optimized up to 300 GeV. The electron 
selection criteria are instead tuned in the multi-100 
GeV range, where the much steeper electron spectrum 
requires an overall hadron rejection power of 1 : 10 4 . 

Events considered for the electron analysis are re- 
quired to fail the ACD vetoes developed to select photon 
events [14]. This removes the vast majority of the po- 
tential gamma-ray contamination. In fact, the geometry 
factor for photons, determined with electron selection 
cuts, is less than 8% of that for electrons at 1 TeV. 
The ratio of photon to electron fluxes is negligible at 
low energies and rises to around 20% at the highest 
energy. This estimate is obtained from a simple, con- 
servative extrapolation to high energies of the EGRET 
all-sky spectrum at GeV energies using a E~ 2A0 power 
law [20]. The overall gamma contamination in the final 
electron sample is therefore always less than 2%. 

EM showers start developing in the LAT TKR, while 
most of the energy is absorbed in the CAL. The mea- 
surement of the lateral shower development is a powerful 
discriminator between more compact EM showers and 
wider hadronic showers. We select on variables that map 
the distribution of TKR clusters around the main track, 
and in the CAL the truncated, second-order moments 
of the energy distribution around the shower axis. A 
further selection derives from the different distributions 
of energy and hits in the ACD between EM and hadron- 
initiated showers. At this stage, the hadron rejection 
power is at the level of 1 to a few 10 2 , improving to 
greater than 10 3 below 100 GeV thanks to the ACD 
selection. 

Similarly to the LAT photon background rejection 
analysis, the remaining necessary boost in the rejection 
power is obtained by combining two probability vari- 
ables that result from training classification trees (CT) 
to distinguish between EM and hadron events [14]. This 
is done using large sets of Monte-Carlo (MC) events 
generated by the accurate LAT simulation package [14], 
based on the Geant4 toolkit [21]. Two CTs are used, 
one built with TKR variables, and a second one based 
on CAL variables, which describe the complete event 
topology. The variables given most weight by the CTs 
are the same or equivalent to those described above. The 
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Fig. 1: Hadron rejection power for the various event 
selection steps. All detector subsystems (ACD, TKR, 
CAL) contribute significantly to the overall rejection. 
The hardening of the rejection power with energy is built 
in the event selection to compensate for the increasingly 
higher proton population at high energies. 



classifiers allow selection of the electrons through a mul- 
titude of parallel paths, each with different selections, 
that map the many different topologies of the signal 
events into a single, continuous probability variable that 
is used to simultaneously handle all valid selections. 
The TKR and CAL electron probabilities are finally 
combined to create an energy-dependent selection that 
identifies electrons with greater efficiency and optimized 
background rejection with respect to a single sequence 
of cuts. The resulting rejection power is flat and better 
than 1 : 10 3 up to 200 GeV and from there rises steadily 
to ~ 1 : 10 4 at 1 TeV in a manner that partially 
compensates for the increasingly larger relative proton 
fluxes with energy, see figure 1 . Conversely, the electron 
selection efficiency, calculated as the ratio of selected 
versus triggered events, has a peak value of 50% at 20 
GeV and steadily decreases down to 12.5% at 1 TeV. 

III. Energy reconstruction and validation 

Energy reconstruction is the other critical aspect of 
this analysis. For EM cascades of several hundreds of 
GeV a large fraction of the energy falls outside of the 
LAT CAL. The shower imaging capability is therefore 
crucial in fitting the longitudinal shower profile in order 
to correct for the energy leakage and estimate the in- 
coming energy with good accuracy. The resulting energy 
resolution for events passing the electron selection is 
shown in figure 2. Since showers are not fully contained 
above 20 GeV, the distribution of the reconstructed 
energy after leakage correction is asymmetric, with a 
longer tail toward lower energies. For this reason we 
quote the full width of the 68% containment of the 
distribution as our energy resolution, and check that 
the full 95% containment does not imply indefinitely 
long tails; see figure 2. Candidate electrons traverse on 
average 12.5 radiation lengths, resulting from the total 
thickness of the TKR and CAL detectors and the effect 
of event selection. 



PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



3 



LAT 95% 

LAT 68% 

q g Q Beam Test MC 68% (60°) 

. A Beam Test Data 68% (60") 

- ....... Beam Test MC 68% (0°) 

< - x Beam Test Data 68% (0 a ) 

£ 0.4- 




earn Test 68% (60") 



10 100 1000 

Energy (GeV) 

Fig. 2: Energy resolution for the LAT after electron 
selection; the full widths of the smallest energy window 
containing the 68% and the 95% of the energy dispersion 
distribution are shown. The comparison with beam test 
data up to 282 GeV and for on-axis and at 60° incidence 
shown in the figure indicates good agreement with the 
resolution estimated from the simulation. 



1000 



800 - 



600 - 



400 




20 25 30 35 40 
Shower transverse size (mr 



50 



Fig. 3: Distribution of the transverse sizes of the showers 
(above 150 GeV) in the CAL at an intermediate stage of 
the selection, where a large contamination from protons 
is still visible. Flight data (black points) and MC (gray 
solid line) show very good agreement; the underlying 
distributions of electron and hadron samples are visible 
in the left (red) and the right (blue) peaks respectively. 



The energy reconstruction algorithm and the event 
analysis rely heavily on the LAT MC simulation. This 
was extensively verified and fine-tuned using beam test 
data for electrons and hadrons up to 282 GeV [22], 
Extensive efforts are made to avoid bias in the event 
selection by systematically comparing flight data and 
MC distributions of likely discriminants of electrons 
and hadrons, and choosing only those that indicate a 
good agreement. Figure 3 shows the very good data- 
MC agreement for the critical variable that maps the 
transverse shower size. 

Systematic uncertainties are determined for all energy 
bins. For each step in the event selection, we scan a 
range of thresholds around the reference value used 
by the cut and derive the corresponding flux versus 
GF curve. We extrapolate the curve to a GF consistent 
with a null cut, and take the relative difference of 
the corresponding flux and the reference as the sys- 
tematic uncertainty associated with the cut. All such 
contributions, taken separately with their signs, and the 
uncertainty of the residual contamination, derived from 
an overall 20% uncertainty in the underlying proton 
spectrum are summed in quadrature. 

The final tuning of the event selection provides a 
maximum systematic error less than 20% at 1 TeV The 
absolute LAT energy scale, at this early stage of the 
mission, is determined with an uncertainty of *^L. 
This estimate is being further constrained using flight 
and beam test data. The associated systematic error is 
not folded into those above as it is a single scaling 
factor over the whole energy range. Its main effect is to 
rigidly shift the spectrum by li,o% w i mout introducing 
significant deformations. 

While event selection is explicitly energy-dependent 
to suppress the larger high-energy background, it is 



not optimized versus the incident angle of incoming 
particles. Nonetheless we have compared the spectra 
from selected restricted angular bins with the final spec- 
trum reported here; they are consistent within systematic 
uncertainties. A further validation of the event selection 
comes from an independent analysis, developed for 
lower-energy electrons, which produces the same results 
when extended up to the the endpoint of its validity at ~ 
100 GeV. Our capability to reconstruct spectral features 
was tested using the LAT simulation and the energy 
response from figure 2. We superimposed a Gaussian 
line signal, centered at 450 ± 50 GeV rms, on a power 
law spectrum with an index of 3.3. This line contains 
a number of excess counts as from the ATIC paper [8], 
rescaled with the LAT GF. We verified that this analysis 
easily detects this feature with high significance (the full 
width of the 68% containment energy resolution of the 
LAT at 450 GeV is 18%). 

IV. Results and discussion 

More than 4M electron events above 20 GeV were 
selected in survey (sky scanning) mode from 4 August 
2008 to 31 January 2009. Energy bins were chosen to 
be the full width of the 68% containment of the energy 
dispersion, evaluated at the bin center. The residual 
hadronic background was estimated from the average 
rate of hadrons that survive electron selection in the 
simulations, and subtracted from the measured rate of 
candidate electrons. The result is corrected for finite 
energy redistribution with an unfolding analysis [23] and 
converted into a flux Je by scaling with the GF. The 
distribution of E 3 x Je is shown in figure 4. 

Fermi data points visually indicate a suggestive devia- 
tion from a flat spectrum. However, if we conservatively 
add point-to-point systematic errors in quadrature with 
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Fig. 4: The Fermi LAT CR electron spectrum (red filled 
circles). Systematic errors are shown by the gray band. 
The two-headed arrow in the top-right corner of the 
figure gives size and direction of the rigid shift of the 
spectrum implied by a shift of of the absolute 

energy, corresponding to the present estimate of the 
uncertainty of the LAT energy scale. Other high-energy 
measurements and a conventional diffusive model [1] are 
shown. 



statistical errors, our data are well fit by a simple 
normalized E~ 3m power law (\ 2 — 9.7, d.o.f. 24). 

For comparison, we show a conventional model [1] 
for the electron spectrum, which is also being used 
as a reference in a related Fermi-LAT paper [24] on 
the Galactic diffuse gamma-ray emission. This uses 
the GALPROP code [4], with propagation parameters 
adjusted to fit a variety of pre-Fermi CR data, including 
electrons. This model has an electron injection spectral 
index of 2.54 above 4 GeV, a diffusion coefficient 
varying with energy as E 1 / 3 , and includes a diffusive 
reacceleration term. As can be clearly seen from the blue 
dashed line in figure 4, this model produces too steep 
a spectrum after propagation to be compatible with the 
Fermi measurement reported here. 

The observation that the spectrum is much harder than 
the conventional one may be explained by assuming a 
harder electron spectrum at the source, which is not ex- 
cluded by other measurements. However, the significant 
flattening of the LAT data above the model predictions 
for E > 70 GeV may also suggest the presence of 
one or more local sources of high energy CR electrons. 
We found that the LAT spectrum can be nicely fit by 
adding an additional component of primary electrons 
and positrons, with injection spectrum J cx tra(E) oc 
E~ la exp{— E/E cut }, E CVLt being the cutoff energy of 
the source spectrum. The main purpose of adding such 
a component is to reconcile theoretical predictions with 
both the Fermi electron data and the Pamela data [7] 
showing an increase in the e + / (e - + e + ) fraction above 
10 GeV. The latter cannot be produced by secondary 
positrons coming from interaction of the Galactic CR 
with the ISM. Such an additional component also pro- 



vides a natural explanation of the steepening of the 
spectrum above 1 TeV indicated by H.E.S.S. data [9]. 
As discussed in [12] and references therein, pulsars 
are the most natural candidates for such sources. Other 
astrophysical interpretations, or dark matter scenarios, 
can not be excluded at the present stage. A more detailed 
discussion of theoretical models can be found in [25], 
as well as in the numerous papers referencing the Fermi 
CR electron spectrum measurement on the arXiv [26]. 
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